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Abstract 



Skew critical problems occur in continuous and discrete nonholonomic 
Lagrangian systems. They are analogues of constrained optimization 
problems, where the objective is differentiated in directions given by an 
apriori distribution, instead of tangent directions to the constraint. We 
show semiglobal existence and uniqueness for nondegenerate skew critical 
problems, and show that the solutions of two skew critical problems have 
the same contact as the problems themselves. Also, we develop some in- 
frastructure that is necessary to compute with contact order geometrically, 
directly on manifolds. 



1 Introduction 

Let M and N be manifolds, suppose /: M — » R is C 1 , and let g: M — > N be a 
C 1 submersion. Given this data, m c G M is a critical point at n G JV if 



This is the standard constrained optimization problem that seeks critical points 
of the objective / subject to the constraint g. 

Appearing in (jl.ip are the derivative of the objective df, the constraint func- 
tion g, and kerTg, which is a distribution on M. Generalizing, we consider the 
data (ct,T>,g), where a is a one- form on M, I? is a distribution on M, and 
g; M — y N is as above. We replace the first condition of with the con- 

dition that a annihilates T>, and we call the result a skew critical problem. 
Skew critical problems occur when an objective function is not differentiated in 
tangent directions to a constraint, but rather is differentiated in the directions 
specified by an apriori given distribution. We are interested in skew critical 



{ 



df(m c )(v) = for all v such that T mc g(v) = 0, 
g{m c ) = n. 



(1.1) 
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problems because, for nonholonomic mechanics, the relevant variational prin- 
ciple is skew [5], and this is also true of the variational discrete analogues of 
nonholonomic systems. 

For mechanics we are interested in existence and uniqueness of skew critical 
problems, by direct perturbation from the point of zero-time change. We have 
a global solution of the (trivial) zero-time problem, and we are interested in 
semiglobal results, which means global along the unperturbed problem, but local 
transverse to that. For discrete nonholonomic systems, we are also interested to 
know that the solutions of two skew critical problems have the same contact as 
the data of the two problems. The skew critical problems of mechanics require 
desingularization at zero-time, essentially by dividing by time. This degrades 
the order matching, which is again recovered by a zero-time symmetry of the 
desingularized problem, and so we must consider the presence of symmetry. We 
are interested in applications to both the continuous and discrete mechanics, so 
we work in an appropriate context of infinite dimensional manifolds. 

In this work, we collect some technical results related to skew critical prob- 
lems. For such a problem (a, V, g), little can be inferred just from the equations 
a(m)\T> — 0, g(rn) — n, without some control imposed on a, T>, and g, so we 
begin in Section [2] with the definition of a nondegenerate skew critical point. 
This corresponds to infinitesimal conditions that, using the implicit function 
theorem, imply there is locally a unique skew critical point for every nearby 
constraint value (Lemma 12. 5[) . If N is paracompact, then a manifold of non- 
degenerate skew critical points along a submanifold Nq C. N can be extended 
along the whole of Nq. We call this result semiglobal because it establishes an 
extension over the whole of Nq, rather that just at one point of N. 

Contact of solutions of skew critical problems is important for discretiza- 
tions of constrained Lagrangian systems, because contact with the exact system 
determines the order of the corresponding numerical methods. Section [3] es- 
tablishes the basic definitions and results about contact. Generally, it often 
happens that cancellations result in one higher contact that would normally 
be expected from data or computation. For example, any Taylor expansion to 
odd order of an even function, is actually the expansion to the next higher or- 
der; a less trivial example is the fact that any odd order self- adjoint one step 
numerical method is one higher (even) order [2]. It is best to understand the 
cancellations geometrically. This kind of "passage to the next order" occurs 
when a geometric object that we call the residual vanishes. In Section[3]we find 
that it is useful to consider the vector bundle analogue of blowing up near the 
zero of a function of a single variable i.e. the function f(t)/t where f(0) = 0. 
The completion of the function is made with the help of the vertical bundle at 
the zero section, and the contact drops by one. We provide, for computing on 
manifolds, Equation (|3.8[) . which computes the residuals of the composition of 
two maps in terms of the residuals of the maps themselves. For skew critical 
problems, it is necessary to consider the contact order of distributions, which are 
subsets rather that maps. This is naturally done using Grassmann manifolds: 
a distribution can be regarded as an assignment of subspaces to base points. 

Finally, in Section |4]wc consider contact for inverse functions and the prob- 
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lems of construction maps from graphs. For graphs, an exchange symmetry of 
the residuals implies that the contact increases by one. In Section [5] we con- 
sider contact for skew critical problems. In the presence of the action of a Lie 
group, we obtain equivariance of the residuals of the skew critical points given 
equivariance of the residuals of the skew critical problems. 

The notations in this work follow those of PQ. We assume without mention 
that the manifolds and submanifolds we use are sufficiently differentiable to 
support whatever operations are involved. 

2 Regular skew critical problems 

Let M and N be Banach manifolds, a be a C k one-form on M, T> be a C k 
distribution on M, and let g : M — > N be a C k submersion i.e. Tg is surjective 
with split kernel. We call (a,V,g) a C k skew critical problem. 

Definition 2.1. A point m c G M is a skew critical point of (a,T>,g) at n £ N 
if 

( a(m c )(v) = for all v G T> mc , 
\ g(m c ) = n. 

A critical point m c of a constrained optimization problem with n = g(m c ) is 
called nondegenerate if the Hessian of is nonsingular. The correspond- 

ing notion for skew critical problems is given below in Definitions 12.21 and 12.41 

Definition 2.2. Let m c be a skew critical point of (a, T>, g). Define the bilinear 
form dx>ci(m c ) : T mr M x T> mc — > R by 

dT>a(m c )(u,v) = (d(ivct)(m c ),u), 

where V is a (local) vector field with values in T> such that V(m c ) = v. The 
skew Hessian of a with respect to g and T> is the bilinear form 

dv, g a(m c ): ker T mc g x P r „ c — > R 

obtained by restriction of dx>a(rn c ). Define dx> , g oi{m c f' : kerT mc g —* T>* m by 

dv, g ct(m c ) b (u) ee d Vtg a{m c )(u, ■). 

Remark 2.3. The definition of dx>a(m c ) does not depend on the extension V: in 
a vector bundle chart of V, the local setup has 

TM = U x (B©F), D = U x (B© {o}), a = a D ©a F , 

where U C B © F is open, a D : [/ -> B* = annF, and a ¥ : U -> F* = annO. 
Supposing that x c G f7 is a skew critical point, two extensions Vi'.U — » B, 
i = 1,2, with Vi(a; c ) = u = V2(x c ) result in iy 1 -v 2 a — (<Xo,Vi — V2). By 



3 



the product rule, the derivative of this at x c is zero since both an and V\ — V% 
vanish at x c , so d{iv 1 ct){x c ) = d(iv 2 a;)(x c ). In contrast to the constrained critical 
problems, skew Hessians are not symmetric since their arguments assume values 
in different vector subspaces. 

Definition 2.4. A skew critical point m c of (a, T>, g) is called nondegenerate if 
d-D,gO.{mcf is a linear isomorphism. 

In finite dimensions, the standard constrained optimization problem (jXTTJ) 
has as many equations for m c as there are unknowns, because g simultaneously 
constrains both v and m c . For the skew problem (|2.I[) . the number of equations 
need not equal the number of unknowns, since g and T> may be unrelated. Defi- 
nition [213] controls this, because if m c is nondegenerate then the fiber dimensions 
of kerT mc (7 and T> mc are equal since kcrT mc g and T > * 1lc are isomorphic. 

Lemma 2.5. Letm c be a nondegenerate skew critical point of a C k skew critical 
problem (a,T>,g), k > 1, and let n c = g(m c ). Then there are neighborhoods 
U 3 m c and V 3 n c such that, for every n £ V there is a unique skew critical 
point m G U of {a,T>,g) such that g(m) = n. Moreover, the map 7: V — > U so 
defined is C k . 

Proof. Using vector bundle charts as in Remark 12.31 the skew critical points x 
such that g(x) = y are obtained by solving F(x) = (0,y), where F(x) = 
(aoix), g{x)) . The derivative of F at a particular x c is 

DF(x c )u = (Da (x c )u, Dg(x c )u) . (2.2) 

The first component is a linear isomorphism on ker Dg{x c ) since x c is nonde- 
generate. Since Dg(x c ) is onto with a split kernel, there is a closed subspace 
K such that D © F = ker Dg(x c ) © K, and Dg(x c )\K is a linear isomorphism. 
From (|2.2p . DF(x c )u = (wi,W2) is continuously inverted by 

u=(Dg(x c )\K)- 1 w 2 , 

u = u+ (Dan(x c )\ ker Dg(x c )) (wi - Da (x c )u), 

and the result follows from the inverse function theorem. □ 

The following semiglobal inverse function theorem is found on page 97 of [4] . 
The semiglobal result for skew critical points which follows that, the proof of 
which is included for completeness, pre-supposes nondegeneracy along a given 
smooth map of skew critical points. 

Theorem 2.6. Let M and N be manifolds and f : M ->■ N be C k , k > 1. 
Suppose that 

1. Mq is a closed submanifold of M , Nq is a closed submanifold of N , and 
/| Mo: Mq — > Nq is a dijjeomorphism; and 

2. f is a local diffeomorphism at every m £ Mq. 
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Then f is a C k diffeomorphism from some open neighborhood U D Mo to some 
open neighborhood V D Nq. 

Theorem 2.7. Let {a,T>,g) be a C k skew critical problem, k > 1, where 
g: M — > N. Suppose that N is paracompact, and that 

a. Mq is a closed submanifold of M, Nq is a closed manifold of N and 
70 : Nq —> Mo is a C k diffeomorphism; and 

b. for all n S No, "fo{ n ) * s a nondegenerate skew critical point of (a,T>,g) 
at n. 

Then there are open neighborhoods U D M and V 3 Nq an d a C k extension 
j:V—>U such that 

1. for all n G V, j(n) is a skew critical point of (a,T>,g) at n; and 

2. 7(71) is the unique skew critical point of (a,T>,g) in U. 

Proof. Applying Lemma [231 at all 70 (no) as n ranges through 7V , there are 
open covers Ui of Mq and Vi of Nq, and C k maps ji : Vi — > Ui, such that, 
for all n £ Vi, Ji(n) is the unique skew critical point of (a, 2?, g) in f/j. By 
shrinking Vi one can arrange 7i(clVi) C Ui where 7, is defined on an open 
superset of Vi. Because is paracompact, its open cover [N \ No, Vi} admits 
a locally finite refinement, so the collection {V,} can be assumed locally finite. 

By Lemma 20.4 of [6], the collection {cl Vi } is also locally finite, so each n G 
[J i V, admits a neighborhood V n that meets only finitely many cl Vi. For each n G 
Ui Vi, the set of indices 

St(n) ee {i : n G c\V % ) 

is finite. No St(n) is empty because every n G lj i Vi is contained in some V and 
hence is in some cl Vi . The set 

V n \ (J{ cl Vi ■■ clVi meets V n and i £ St(n) } (2.3) 

an open neighborhood of n because it subtracts from V n only finite many closed 
sets, and it has the property that if any of its members is in any cl Vi then 
i G St(n). Replacing each V n with (|2.3|) . it can be assumed that St(n') C St(n) 
for all n' G V n . 
Defining 

(J (f'Wn fl u *)> V = \J 7 r\u), 

«eUj Vi i£St(n) ' i 

we can show the following facts. 

1. Mo C U: if m G Mq and n = g{m) then m G g _1 (n) and n G clVi for all 
i G St(n) so 7,(n) G ^(clVi) C [/, for all i G St(n), hence m G t/. 



5 



2. U is an open neighborhood of M : if m £ U and n = g(m) then 



me<rV)n f) c/.cf'fKln f| t/i. 

iGSt(n) ieSt(n) 

The last set is open because it is the intersection of finitely many open 
sets. Also, 



g-\v n )n fl u i= |J (ff-V)n f| ui) 



n'eV n v ieSt(n') 

c u. 

Thus there is an open neighborhood of m that is contained in U. 

3. U has the property that, for all mi, m2 G £/, g(m{) — g(jri2) implies that 
there is an i such that mi and m2 are both in ?7j. Indeed, any such mi 
and m2 are members of 

. 9 -V)n n u *> 

ieSt(n) 

where n = g(m\) = g(m 2 ), and so both mi and m 2 are members of any 
Ui for any i £ St(n). 

Let n e y. Then n £ "f^ 1 ^) for some z and m = 7,(n) is a skew critical 
point of (a, T>, g) in J7. If m' £ U is another such skew critical point then 
g(m) = g(m'), and m and m! both lie in a single Uj. By definition of the Uj 
there is only one skew critical point of (a, V, g) in Uj, som = m'. Thus for all 
n £ V there is a unique skew critical point of (a,T>,g) in U. Define 7: V — > U 
by this correspondence. By the uniqueness used to define 7, the restriction of 
7 to any 7 2 rl (t/) is 7$, which is C k , and the 7 J ~ 1 (£7) cover V, so 7 is C k . □ 



3 Order Notation and Residuals 

Given two functions fi(x), i = 1,2, of a single variable x £ M, the standard 
definition of fi(x) = f2(x) + 0(x r ) is that there are numbers 5 > and C > 
such that |/i(x) — /2(^)| < C|a:| r for |x| < 5. If the functions /, are C r , r > 1, 
then /i(a;) = /2(a;) + 0(x r ) if and only if there is a continuous function, say 
6f(x), such that fi(x) = /2(x) + x r Sf(x). The following definitions export the 
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second formulation to the context of manifolds. 
Definition 3.1. 

1. Let M be a manifold and Hm ■ M — > K be a C°° function which has as 
a regular value. The pair (M, Hm) will be called a manifold. 

2. Let (M,Iim) and (JV, /ijv) be manifolds. A C fc mapping f: (M, ftjif) — > 
(JV, /ijv) is a C fc mapping / : M —> N such that /ijv ° / = ZiAf- 

3. A mapping /: (M,}im) — > N or f : M — > (AT, ft,jy) is a mapping 
f:M^N without any conditions involving /im or /ijy. 

Definition 3.2. Let (M, %/) and iV be manifolds and /, : (M,Hm) — ► AT, 
z = 1,2, be such that /i = /2 on /iT, (0). Define /2 = /i + 0(h r M ), r > 1 if, for 
all mo S /i^ (0), there is a chart ^ at rto = /i(mo) € AT, and there is a function 
(<5/)„ defined near mo, and continuous at mo, such that 

f (/a( m )) - v(fi{m)) = h M {m) r {5f) y {m), 

for all m in some neighborhood of mo- 

The definition of /2 = /i + 0(h r M ) docs not depend on the coordinate 
chart: if and £ are two coordinate charts at no, as in Definition 13 - 21 and 
for m near to mo, 

v{h{m))-v(h{m)) 

= {vo V - x )(y{h{m))) - (PoOK/iM)) 

= (P o z/" 1 ) (i/(/i(m)) + h M {m) r {5f) v {mj) - (i> o i/" 1 ) (i/(/i(m))) 



o cft v 
h M (m) 



Do V - 1 )L(f 1 (m))+th M (my(Sf) u (m)) dt 



D^ov- 1 ) Ufi(m))+th M (m) r {5f) v {m) ) dt 



as required. 

The quantities (Sf) u (mo) and (Sf)o(n%o) transform as tangent vectors. In- 
deed, 

h M {m) r {5f)v(jn) = vfaim)) - v(fi(m)), 



so 



D{v o v- x ) Uh{m)) + th M (m) r (Sf) u (m)) dt 



At mo G h^j (0), and setting no = fi{n%o), 

(Sf)z(mo) = D(Pov- 1 )(v(n )){Sf)„(mo), 

as required. 
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Definition 3.3. Let (M, Jim) be a manifold, f 2 = fi+0(h 7 M ), and m € h^(0). 
The vector res r (f 2 , /i)(m) £ T n N with representation (8f) u (m) for any chart ^ 
is called the r-residual of f 2 with respect to fi . 

The residual res r (/2,/i) is defined only on h~^~(0) C M and takes values 
in TN . The condition f 2 — fi + 0(h r M ), can be localized to a point of M 
or a subset of M in the obvious way, and the residual will be correspondingly 
localized. In general, jets of mappings between manifolds carry an affine action 
by a geometrically based vector space, amounting basically to the first nonzero 
term of the Taylor series of the difference between two mappings. Also, the 
notion of contact below is the same as the contact equivalence in the definition 
of jets [3]. 

If (M, Iim) is a manifold, then, since is a regular value of Hm, there are 
/im -adapted charts at each m £ ^m(O) i.e. charts such that the local rep- 
resentative of h\i is the projection (x,t) i— > t. We can prove an equality or 
formula concerning residuals in any chart since residuals are geometric, and in 
particular, we can always use an /iM-adapted chart. 

Suppose that fi: (U,h(j) — ► V C F, i = 1, 2, are C r , where U is an open 
subset of E x R, where E and F are a Banach spaces, and where hjj{x,t) = t. 
For fixed x, the Taylor expansions in t about t = of the fi are 



8 f t r 8 r f- 

fi(x, t) = fi{x, 0) + t-jg(x, 0) + • • • + - { -Q]r(x, 0) + i^Or, t) f 



where [T] 



R rti (x,t) = j ( r _i)i ^atr^'**)" aF (a: ' 0) ) ds - 



The condition that fs = fi + 0(h\j) at (x, 0) is thus equivalent to the condition 
that these Taylor expansions match at (x, 0) up to and including the degree r — 1 
term. So, given this, 



h{x,t) - h(x,t) 

t r ( d r f 2 



(X, 0) - ^(X, 0)) + Rr, 2 (x, t) t r - Rr A (x, t) f , 



r! \ 9r 

which identifies (5/)„(x, t) in these coordinates as 

(6fUx, t) = i (^(», 0) - ^(z, 0)) + fl^x, t) - iZr.iCa:, t). 
Setting t = 0, the residual is 

(/ 2 (M)-/iM))- 



If res r (f 2 , /i) = 0, then the Taylor series of /i and / 2 agree up to and including 
terms of degree r, one more than the degree r — 1 agreement implied by f 2 = 
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fi + 0{Kjj). Thus, if it is necessary to establish with some computation, that 
two functions which differ at order r, actually differ at order r + 1, then one can 
accomplish this by showing that res r (/2, /i) = 0. 

Definition 3.4. If M, h,M, and fi are as in Definition 13.21 then f\ and fi have 
order ft77 contact, or just have contact r — 1, if fa = f\ + 0(h r M ). 

If 7r : E — > M is a vector bundle then kerT7r is a subbundle of TE, called the 
vertical subbundle. Recall that the tangent space of the zero section defines a 
natural horizontal subspace, so any vector of TE at the zero section can be split 
into horizontal and vertical parts. This splitting can be defined by the linear 
isomorphism 



TM © .E — > TE : (v m , w m ) ^ ^~ 

at 







d 



'm(t) 



tw m , (3.1) 



where m(t) is a curve in M such that m'(0) — v m . If z £ To m E then denote the 
horizontal and vertical parts of z by horz £ T m M and vert z £ E m , respectively 
i.e. the inverse of (|3. 1|) is z t— > (hor z, vert z). 

If / is a C 1 function such that /(0) = 0, then it is elementary that 

MO 
l/'(0), < = o, 



is continuous. The purpose of Lemma 13.51 is to show that a mapping on a 
manifold can be smoothly divided by a real function that takes values in a 
vector bundle and is in the zero section 0(E) if the function vanishes. 

Proposition 3.5. Let (M,Hm) and N be a manifolds, and let it: E — » N be 
a vector bundle. Suppose that f:U — > E is C k , k > 1, and t/iat /(m) G O(-E) 
whenever = 0. T/ien /or a/Z to swc/i f/iai hM(m) — 0, f/iere is a unique 

e(m) £ E n ftf m \) such that 

vertT m f(v m ) = (dh M (m)v m )e(m), v m £ T m M. (3.2) 

Moreover, the function f : M — > E defined by 

/(to) = < h M (m) 

[e(m), h M {m) = 0, 

is C k - 1 . 

Proof. If /ijvf(m) = and u m £ kerd/ijv/(m) then u m = c'(0) for some curve 
c(t) € hjj(0). Since / is in the zero section wherever %m is zero, it follows 
that / o c[t) takes values in the zero section, so (/ o /im)'(0) is horizontal. Thus 
vertT m /(w m ) = for all v m £ kerdft,j\/(m), so there is a unique e(m) 6 E^(tr m )) 
satisfying (|3.2p . 
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We can set up an /ij^f-adapted chart { (x, t)} on M and a vector bundle chart 
on N, so that E = {(y,e)}, and f(x,t) = (fo(x,t), fi(x,t)). Then /i(x,0) = 
for all x, so 

veTtT m f(x,0)(Sx,6t) = ^(x, 0)Sx + ^(x, 0)St = ^M)« ( 3 - 3 ) 

and 

dh M (x,t)(Sx,St)=6t. (3.4) 
By comparison of (|3~2f with (|3~3|l and ([3~4]) . 

e(x,0) = ^(x,0), 



and it is required to show that fx defined by 

fi(x,t) 

fi(x,t) 



t 

df: 



t^O, 



T^M), * = o, 

is C k ~ 1 . At any (xo, 0) the Taylor expansion of /i is 

fi(x,t) = Df 1 (x o ,0){Sx,t) H 

h ^ fe /i(^o, 0)(<5x, i) fc + i? fc (x, i)(fa, i)* 



(3.5) 



where 5a; = a; — xo and R(xq,0) — 0. By differentiating /i(x, 0) = in i, 
L> l /i(^o,0)((5a;,0) i = for 1 < i < fc, and substituting t = into gives 
i?(x, 0)(5x, 0) fc = 0. Thus the left side of (|3.5|) has t as a factor and 

/ifot) = ^(ar o ,0) + -D 2 h(x Q ,Q){5x : t) 2 + ■■■ 

(3.6) 

• ■ ■ + ^> fc /i(zo,0)(t,/i)* + ^(z, *)(&,*)*. 

Each of the first k terms of (|3.6p are polynomial in (5x,i) and the remainder 
is polynomial in [5x, t) of degree k — 1 with coefficients functions of (x, i) that 
vanish at (xo, 0). Thus, by the converse of Taylor's theorem p], /i(x, t) is C^" 1 
at any (xo , 0) . □ 

Proposition 3.6. Let (M, Km ) a^rf N be a manifolds, let ir: E ^ N a vector 
bundle, and suppose fi and fi are as in Proposition \3.5l with k > r. Then 
h = A + O^m 1 ) ifh = fi + 0(h r M ), r > 2.^ Moreover, res r (/ 2 ,/i) tafces 
values in the vertical bundle of E and res r_1 (/2, /i) = res r (/2,/i)- 
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Proof. Assume the context and notations of the proof of Proposition ^ .51 so that 

fl(x,t) = (fi,o(x,t),fx,i(x,t)), fs(x,t) = (f2,o(x,t)J 2 ,i(x,t)). 

Since f 2 ,o(x,t) = f± t o(x,t) + 0(t r ), the r — 1 residual of the first components 
of f% and f\ is zero, and it suffices show that 

f 2 ,i(x,t) = ft^x^) +t r - 1 6f(x,t) (3.7) 

given / 2i i = fi t i + t r 5f, where Sf is continuous and 



fi,i(x,t) = 



t 

dfi, 



-(x,0), t = 0. 



Of 

Equation (|3.7|) can be shown in the two cases t = and t ^ 0: For < = 0, 
— (z,0) - -^-M) = lrm— — = hmt Sf(x,t) = 0, 

so even / 2j i(cc, = fi,i{x, t) in this case, whereas for t ^ 0, 

? / ,x /2,1PM) f hl (x,t)+t r Sf{x,t) f .s 
J2,i{x,t) = = = Ji,i(a;,tJ + t of(x,t). 



□ 



Proposition 13.71 is a key result because it can be used to compute residuals 
without the invocation of local charts. Note that if (M, 1t,m) and (N,Hn) are 
manifolds and /: M — > TV is a C 1 function such that f(h~^(0)) C /i^ (0), then 
for all m £ h~^(0), d(hjv f)(m)v m = for all u m such that dhM(ni)v m = 0. 
So one can define /: /^/(O) K by 

f)(m) = f{m) dh M {m). 



This is an instance of Proposition 13.51 and it follows that h^j defined by ex- 
tending d(hpf o f)/fiM to / on hpj(0) is continuous. 

Proposition 3.7. Let (M, Km), (N,hM), and P be manifolds, and suppose 
fi : M -> N and g { : N -»■ P, i = 1, 2 are C 1 and satisfy ^(h^iO)) C /^(O), 
/2 = /1 + 0(h r M ), and g 2 = gi + 0(h r N ). Then g 2 o / 2 = 51 o / x + 0{h r M ). 
Moreover, if '/(^(m) = cmd n = fi(m), then 

res r (g 2 ° /2,5i /i)M = /2(m) r res r (g 2 , 5i)W + T n g x res r (/ 2 , fi)(m). (3.8) 
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Proof. It suffices to consider the local setup where E, F, and G are Banach 
spaces, [/CE and V C F are open, / : U -> V, .9 : V —> W, hjj : U — > K, and 
/ly : V -> E. Then 

(.92 o / 2 )(x) = gi(fi(x) + hu(x) r 5f(x)) + h v (f 2 (x)) r Sg(f 2 (x)) 



Dg^hixJ + shuixYfix^ds 



ds 

+ h v (f 2 (x)) r 5g(f 2 (x)) 

= gi(h(x)) + hu{x) 

+ h u (x) r h V j 2 (x) r 5g(f 2 (x)). 
Assuming x satisfies hjj(x) — and putting y = fi(x), results in 

res r (g 2 0/2,510 h)(x) = D gi (y)Sf(x) + h VJ (x) r 5g(y), 
which is the local form of 



Sf(x) 



□ 

Remark 3.8. If r > 2, then T n gi and f 2 can be replaced by T n g 2 and /1 respec- 
tively in Equation (13. 8|) . If h,N°fi — hu then /j = 1 in any case. Also, if g\ = g 2 , 
one can dispense with fiN and the assumption that fi(hjl(0)) C ft,^- (0), ob- 
taining the formula 

res r {g h, 9 f\){™) = T n g res r (f 2 ,fi)(m). 

Let f : £ -4 B be a vector bundle with typical fiber E, and let Eo be a 
closed split subspace of E. We will have use of the Mo-Grassmann bundle of E, 
denoted 7r^ Eo ' S ' : G(E , E) — > M, by which we mean the set of subspaces of the 
fibers of E that are linearly isomorphic to Eo i.e. the coset space of continuous 
linear injections (with closed split image) of Eo into the fibers of E by the action 
of GL(Eq). The projection nfj- E °' E ^ associates subspaces of the fiber E m to m, 
and the typical fiber of G(Mq,E) is the Grassmann manifold Ge (E). For more 
information on Grassmann manifolds in the Banach space context, see pQ. 

Remark 3.9. If E' is a C r subbundle of E, with typical fiber E', then there 
is defined the C r map Le>'- M — > G(W,E) that assigns to any m S M the 
subspace L^'(m) — E' m . 

Remark 3.10. As is well known, the tangent space at B G Ge (E) is canonically 
hom(B,E/B). Indeed, if M(t) is a C 1 curve in G Eo (E) with B(0) = B, then 
choose a splitting E = B © F and define B: B -> F by 



d_ 

di 



where 7r^/ B denotes the projection of E to the quotient E/B, and 7rJ® F denotes 
the projection to B using the decomposition E = B © F. One verifies that B is 
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independent of the choice of the complement F and depends only on the tangent 
vector of B(t) at t = 0. 

Definition 3.11. Let {M, Km) be a manifold and let T>i, i = 1,2 be distri- 
butions on M such that V]\h^}{Q) = V 2 \h^(0). Define V 2 = V 1 + 0(h r M ) if 
iB 2 = tDi +0(h r M ), and define res r (X> 2 , X>i ) = res r (t D2 , i Dl ). 

In the context of Definition 13.111 note that the assignment of the fibers of 
subbundles into the Grassmann bundle preserves fibers, so 7T^" oi Dj = 1 A /, 
and 

rp G„(TM) ( \ / G„(TM) G V (TM) \ 

Ttt m res(t C2 , tuj = res(?r M 'oi D2 , r M o i23l j 
= res(l M , 1m) 
-0, 

which shows that, for all m 6 hj}(Q), res r (T>2, T>i)(m) is a vertical vector in 
T (G-d(T M)j . Such vertical vectors are derivatives of curves in the corresponding 
fiber i.e. derivatives of curves in the Grassmann manifold Go(T m M). Thus, the 
residual of two vector bundles T>i at m is an element of the tangent space at 
the common element (2?i) m of the Grassmann manifold Gn{TM), which, by 
Remark 13.101 can be regarded as an element of hom{(T>i) m ,T m M/(T)i) m ). 

4 Equations 

Computing with the order notation on manifolds might require the determining 
the contact or residual of the solutions of two implicit equations with a given 
contact or residual. A most basic result that enables this sort of argument is 
Proposition 14.11 which guarantees the contact of two inverse mappings, given 
the contact of two diffeomorphisms. 

Proposition 4.1. Let (M,Hm) md (N,Jin) be manifolds, and let /j: M N 

be C k diffeomorphisms, k > 1, be such that fi maps h~^j(0) into hjj 1 ^), i = 1,2. 
Then h = fx+ 0(h r M ) implies f^ 1 = /f 1 + 0(h r N ). 

Proof. Suppose / is such that g 2 ° f% = gi ° fi + 0(h l M ). This is true for I = 1, 
because g 2 o / 2 = g% o f 1 on h]^(0). Taking the residuals of fi 0^ = 1, one 
obtains 

= rcs'(/ 2 °92, fx °9i)(n) = g 2 (n) 1 res'(/ 2 , /i)(m) + T n fi res l (g 2 ,gi)(n). 

Thus ves l {g 2 , gi )(n) = if res'(/ 2) f x ){m) = i.e. g 2 = 9l + 0(h l N ) if f 2 = 
fi + 0{h l M ), which inductively gives g 2 = g\ + 0{h r N ). □ 

Another requirement is to semiglobally construct mappings from graphs. 
Proposition 14.21 uses the semiglobal inverse function theorem to provide such a 
result for a perturbation of an identity mapping. 
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Proposition 4.2. Let M and (N, h N ) be manifolds. Let 7: U C N M X M 

be C k , k > 1. Suppose that /i^(0) C 17 anc? 7 1 ft,jy (0) is a diffeomorphism 
to A(M x M). Tften there are neighborhoods U C {7 o/Zi^O) and V C M x K 
0/ M x {0} smc/i fftai, /or aZZ (m,h) £ V, there is a unique rh € M sucft i/iai, 
/or some n EU, 7(71) = (to, to) and /i/vf/i) = h. The map / 7 : V M defined 
by f-y(m, h) = rh is C k . 

Proof. Let 7Ti and 7T2 be the projections on M x M i.e. 7Tj(mi,m2) = m„ 
i = 1,2. Define ip: U -> M x M by t/>(?7) ee ((tti o^n),^ (n)). The map 
t/i is a diffeomorphism from ^^(O) to M x {o} and, by the inverse function 
theorem, is a local diffeomorphism at each point of hj^fi). By Lemma |2~51 if) is 
a diffeomorphism from a neighborhood U C IT of /i^- 1 (0) to a neighborhood V 
ofMx{0}. 

If (m,h) E V, then let n S U be such that ^>(n) = i m >h), and define 
m = 772(7(71)), so that / 7 (m, h) = rh = (n20'joip~ 1 )(m, h). From ^(n) = (m, /i) 
follows (711 (7(77)) , /ijv(n)) = (m,h) so 7(71) = (to, to) and h^{n) = ft, which 
are the required properties of to. If there is another such, say rh', then there 
would have to be an n' 6 U such that 7(77') = (m,fh!) and hjv(n') = h, so 
rjj(n') = (m,h) = ip(n) which, since ip is a diffeomorphism, implies rt = 77'. 
Hence (777,777) = 7(77) = 7(71') = (to, to'), so 777, = 777'. □ 

Proposition 14.21 establishes the contact of the mappings constructed from 
graphs is equal to the contact of the graphs. Further, the mappings have one 
higher contact if there is present a symmetry condition for the residuals of the 
graphs. 

Proposition 4.3. Let [M, Hm) and (N, Hn) be manifolds and 7* and fi be as in 
Proposition^- S\ Then / 72 = / 7l +0(h r M ) if 72 = 71 +0(h r N ). If res r (j2, 71) is 
symmetric i.e. 8j (n) = S~f 2 (n) for all n € N, where res r (72,7i) = (07 1 , S~f 2 ), 
then / 72 = / 7l +0{h r+1 ). 

Proof. Assume the context and notations of the proof of Proposition ^. 21 Since 
fji = 7T3 o 7j o if)7~ , where ^ = (72, h^), Propositions 13.71 and 14. II imply / 72 = 
/ 71 + (9(/7^) if 72 - 71 + 0{h r ). Then 

so, taking the residuals of this equation at 77 S ft.^ (0), and setting to ee 
^lfji («•)): gives 

r( m , m )7r 2 res r (72,7i)(77) 

= res r (/ 72 o (tti o 72, /l/y), / 71 o (7Ti o 7l , hjv)) (77) 

= res r (/ 72 , / 7i )(tt7, 0) + T (mi0) / 7l res r ((7Ti o 72, /ijy), (tti li,h N )) (n) 

= res r (/ 72 , / 7i )(t77, 0) + T (m:0 )/ 7l (Tttx res r (7 2 , 71), 0)(n). 
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Also, / 7l (tti (mi, 7712), 0) = uii for all mi 6 M, so the last term of the equation 
immediately above is T( m m )7ri (res r (72, 7i)(«)) , and hence 

res r (/ 72 , / 7l ) (to, 0) = T (mim) 7r 2 res r (72 , 71 ) (rc) - T (m , m) 7ri res r (72 , 71 ) (n) , 

which is zero if res r (72, 7i)(n) is symmetric. □ 

5 Skew critical problems 

Theorem l5.1l is a main objective of this work. It uses the infrastructure we have 
developed to show that the contact of solutions of nondegenerate skew critical 
problems is the same as the contact of their data. Moreover, the residuals of 
the solutions are determined geometrically through the residuals of the data. 

Theorem 5.1. Let (M, Jim) and (N^h^) be manifolds and suppose a 1 , gi, 7, 
andVi, i — 1,2 are as in Theorem \2. 7| and Mq C h^j(0), No C /^(O). If a 2 — 
a l +0{h r M ), 92=gi + 0(h r M ), andV 2 =V 1 +0(h r M ), then 72 = 71 + 0(h r N ). 

Proof. It suffices to consider the local setup at x — x c , where 

1. x c € U C E and VCF are open in Banach spaces E and E, respectively, 
and g(x c ) = y c ; 

2. the fiber of T>i at x is the graph {e + Aj(x)e : e e D}, where E = BeB- 1 , 
Ai(x) : B — > B 1 - , and Ai(O) = 0; 

3. a': {/-> E*. 

In this setup, x = Ji(y) are determined by the equations Fi(x) = (0,y) such 
that Fj : [/ — > B* x F is defined by 

Fi(x) = (a Ai (x),g 4 (x)), QfA; = ct{x) o (t DE + Ai(x)), 

where tjjE is the inclusion of B into E. The domain of the Fi has the lo- 
cal representative hjj of Km, and the codomain of the Fi has the function 
^D*xv(a,y) = hy{y) where hy locally represents h^. Fi = F\ + 0(hy) since 
a 1 = a 2 + 0(h[j) and Ai = A 2 + O(h^), and since composition of linear maps 
is continuous and bilinear. Also, Xq IS cL nondegenerate skew critical point for 
both problems corresponding to i = 1, 2, so each Fi is a local diffeomorphism 
at x c - From Proposition 14. 1[ and near (0, y c ), F^ 1 — F^ 1 + 0(h^ s » xV ), which 
from 7,(y) = (0, y) implies 71 = 72 + 0(/iy), as required. □ 

In the context of Theorem 1 5. 1[ we will need to know that the residuals of the 
solutions 7j depend only on the residuals of a 1 , T>i, and gi. For this, it suffices to 
show that, given a skew critical point m c £ M c at n c £ No, u c = res r (72, 7i)(n c ) 
is the unique solution of F(u) = subject to the constraint G{u) = 0, where 
F: T m M -> V* ma is defined by 

F(u) = 7 2 (n c )(iDa :L (m c )(M) 

(5-1) 

+ res'" (a , cr)(m c ) + a (m c ) o res r (2?2, £>i)(to c ) 
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and G : T m< M -» T n N is denned by 

G(u) = Tgi(m c )u + res r (g 2 , gi){m c ). (5.2) 

Here a 1 (m c ) annihilates (£>i) mc and so descends to a (m c ) in the quotient 
T mc M/(2?i) mc . If r > 2 then the index 1 occurring asymmetrically in (|5.ip 
and (|5.2p , such as in the fragment d-oa 1 {m c ) , can be replaced by the index 2 
because the data of the skew critical problems are assumed to match to order. 
To show (|5.1I) and (|5.2p . note that, in the local setup, a vector field in T>i 
extending any e £ D is available as x i— » e + Ai(x)e, and so 



d Dl a 4 (x c )(u, e) = 4" 



(a 8 (x c + ut), e + Ai(x c + ut)e) = (DctAi (x c )u, e). 



t=o 

Since a Ai (7i(y)) = 0, the residuals of this for i = 1, 2 are zero, so 

= 72(yc)res r (aA 2 ,aA 1 )(£ c ) + Da Al (x c ) res r (7 2 , li){y c ) 
= 72(2/ c )res r (a 2 ,a 1 )(x c ) +a 1 (x c ) ° res r (A 2 , Ai)(x c ) 
+ dr> 1 a 1 (a; c ) b res r (72, 7i)(2/c), 



which is the local version of (15. ip . The constraint (15. 2p follows from the equa- 
tion rcs r (.92 o 72,5i ° 7i)(«c) = 0, since #i oj^n) = g 2 o j 2 (n) = n. 

Suppose one has skew critical problems as in Theorem l5.H where the unper- 
turbed problem is equivariant under the action of a Lie group. For the applica- 
tion we have in mind, Q is not a symmetry group of the full critical problem: 
only the residuals of the unperturbed problem are equivariant. Then Proposi- 
tion [572] below shows that the residuals of the solutions are equivariant. Recall 
that, if a Lie group acts on a manifold M, then it acts by lifts on TM and T*M, 
and also in the obvious way on the vertical bundles of TM and T*M, and on 
any Grassmann bundle of TM. 

Proposition 5.2. Let (M,Hm) and (iV, hrj) be manifolds suppose a 1 , gi, ji, 
and T>i, i = 1,2 are as in Theorem \5.1i Suppose that a Lie group Q acts on M 
and N , and 

1. g t : (M, h M ) -> (N, Iin) i.e. h N o g, t = Km; 

2. T>i\Mq are tangent to h~^(Q) and are invariant, a l |T(/V^- (0)) are invari- 
ant, and gi\h~7j (0) are equivariant; 

3. res 1 * (a 2 , a 1 ), res'* (52, 9i)> oxid res r (72, 71) are equivariant. 

Then Km °7i = hjsi, and res(72,7i): (0) — » T(/i^ (0)) is equivariant. 
Proof. Since gi o 7i(n) = n, 

h N (n) = h N (gi o %(n)) = h M ° Ji(n) 
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and Km °1% = hjsr follows. Fix r € Q and let a 1 = r*a, T>i = rT>i, and cji = r* gi, 
where r* denotes pull-back by m t— » rm. Note that the maps % = Tji give the 
skew critical points of (a 1 , X>i, <?i). By (|5.ip and (|5.2p , the residuals res r (72,7i) 
and res r (72 , 71 ) are determined by the residuals of the data in the skew problems 
(a 1 , T>i,gi) and {a 1 ,T>i, cji), respectively. So the assumed equivariance of the data 
residuals implies that the residuals res r (72,72) and res r (72,7i) are equal, and 

rres r (72,7i)(n) = res r (r7 2 , Tji)(rn) = res r (7 2 , 7i)(rn) = res r (7 2 , 71) (to), 

as required. □ 
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